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Abstract 

The  steady-state  and  dynamic  performance  of  Hybrid  ARQ  protocols  in  slot¬ 
ted  direct-sequence  code  division  multiple-access  networks  is  analysed.  A  Markov 
model  is  used  to  derive  throughput-delay  expressions  in  terms  of  the  channel  cutoff 
rate  and  capacity.  Network  design  parameters  are  identified  and  their  dependency 
on  systems  parameters  is  examined  in  det^.  It  is  shown  that,  for  a  given  number 
of  users,  traffic  intensity,  and  bit  energy-to-background  noise  ratio,  there  is  an  op¬ 
timal  probability  of  retransmission,  code  rate,  and  processing  gain  that  maximises 
performance. 


I  Introduction 


Code  division  multiple-access  (CDMA)  continues  to  receive  interest  as  a  viable  network 
multiple  accessing  technique  for  both  civilian  and  military  applications.  This  interest  is 
prompted  by  the  numerous  advantages  which  are  afforded  by  its  use,  including  decentralized 
control  for  uncoordinated  access,  zero  access  delay,  gradual  degradation  in  performance  with 
increased  number  of  users,  increased  noise  immunity,  message  privacy,  and  anti-jam  (AJ) 
and  low-probability  of  intercept  (LPI)  protection.  While  most  of  the  work  to  date  has 
dealt  with  the  physical  level  issues,  a  number  of  papers  which  treat  the  issues  at  the  link 
and  network  levels  have  recently  appeared  in  the  literature.  In  [1],  Raychaudhuri  analyzed 
the  CDMA  link  performance  of  a  paired-off  user  topology  for  several  packet  generation 
models.  More  recently,  Polydoros  and  Silvester  [2]  developed  an  analytic  framework  in 
order  to  identify  those  design  parameters  which  summarize  the  effect  of  various  network 
characteristics  on  link  level  performance.  In  a  paper  related  to  the  present  paper,  Stiiber 
[3]  developed  an  analytic  model  for  evaluating  the  performance  of  slotted  ALOHA  direct- 
sequence  CDMA  networks  with  forward  error  correction  (FEC).  In  this  paper,  an  approach 
similar  to  [3]  u  followed  to  investigate  the  link  level  issues  associated  with  the  use  of  Hybrid 
Automatic-Repeat-Request  (ARQ)  protocols  in  slotted  direct-sequence  CDMA  networks. 

The  key  performance  measures  of  random  access  networks  are  their  throughput  and 
delay  characteristics.  Hybrid  ARQ  schemes  are  often  employed  to  improve  these  character¬ 
istics  in  terms  of  increased  throughput  and  network  reliability.  The  purpose  of  this  paper  is 
to  develop  an  analytical  framework  for  evaluating  and  optimizing  the  throughput-delay  per¬ 
formance  of  Hybrid  ARQ  protocols  in  slotted  direct-sequence  CDMA  networks.  Asymptotic 
performance  bounds  are  derived  in  terms  of  the  channel  cutoff  rate  and  capacity.  These 
bounds  are  considered  asymptotic  because  they  assume  the  use  of  error  detecting/correcting 
codes  operating  arbitrarily  close  to  the  channel  cutoff  rate  or  capacity. 

An  overview  of  this  pap>er  is  as  follows.  In  section  H,  the  network  model  is  defined 
in  terms  of  its  physical  and  link  level  characteristics.  A  Markovian  model  is  employed 
in  section  HI  to  derive  asymptotic  throughput-delay  expressions  in  terms  of  the  channel 
cutoff  rate  and  capacity.  In  section  IV,  numerical  results  are  presented  to  demonstrate  the 
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asymptotic  performance  of  the  CDMA  network.  Extensive  analysis  is  made  to  investigate 
the  dependency  of  the  probability  of  retransmission,  code  rate,  and  processing  gain  on 
system  parameters  such  ^as  the  traffic  intensity,  number  of  transceivers,  and  the  level  of 
background  noise.  Although  the  effects  of  jammer  noise  are  not  considered  here,  they 
provide  the  motivation  for  this  paper  and  are  addressed  in  another  article  [4].  The  results 
of  this  paper  form  the  ‘baseline  pei-foriucnce’  for  the  ismmiug  analysis  in  [4].  Operational 
considerations  are  presented  in  section  V.  In  the  appendix,  the  state  transition  probabilities 
are  derived  for  the  Markov  Model  of  section  m. 

II  Network  Model 

The  network  consists  of  U  radio  units  arranged  in  a  paired-off  user  topology  similar  to 
that  discussed  in  [1,2].  With  this  topology,  (7/2  =  Nt  =  Njt,  where  Nt  is  the  maximum 
possible  number  of  active  transmitters,  and  Nji  is  the  maximum  possible  number  of  active 
receivers.  To  simplify  the  notation,  the  population  size  will  be  expressed  in  terms  of  N,  the 
maximum  possible  number  of  active  receiver-transmitter  pairs  during  a  given  slot  period, 
where  N  =  Nt  =  Nr,  With  this  topology,  only  one  packet  may  be  directed  to  a  given 
receiver  during  a  channel  slot  and  there  is  complete  receiver  availability  whenever  the  cor¬ 
responding  transmitter  is  active.  Therefore,  there  is  no  topological  competition,  and  the 
channel  is  not  a  capture  channel.  The  only  disruptive  interaction  .between  users 

is  multiple-access  interference. 

Radio  units  access  the  network  by  using  direct-sequence  CDMA  with  a  standard 
slotted  ALOHA  protocol.  Information  is  transmitted  in  the  form  of  packets,  one  packet  per 
time  slot.  Packet  flow  for  the  network  is  shown  in  Fig.  1.  Each  of  the  network’s  transceivers 
can  be  in  one  of  two  modes;  origination  or  blocked.  In  the  origination  mode,  the  probability 
of  transmitting  a  packet  in  the  time  slot  is  geometrically  distributed  with  parameter  po, 
where  po  is  the  probability  that  the  user  will  transmit  the  packet  in  the  next  time  slot. 

Prob  {  packet  transmission  in  the  time  slot  }  =  po(l  -  Po)*“^  (1) 

When  either  multiple-access  interference  or  background  noise  causes  a  packet  to  be  received 
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in  error,  a  user  enters  the  blocked  mode.  In  this  mode,  the  probability  of  retransmission 
is  also  geometrically  distributed  but  with  parameter  Pn  where  pr  is  the  probability  that  a 
user  will  retransmit  its  packet  in  the  next  time  slot.  While  a  user  is  waiting  to  retransmit, 
it  is  considered  blocked  or  backlogged  because  it  cannot  transmit  (but  may  receive)  a  new 
packet  until  the  retransmitted  packet  is  received  without  error.  In  practice,  limits  are  set 
on  the  maximum  number  of  retransmissions  allowed  for  a  given  packet.  Here,  we  define 
the  traffic  intensity  as  the  average  number  of  packets  transmitted  per  slot,  and  denote  the 
new  and  retransmitted  traffic  intensities  as  Vo  respectively.  At  the  input  to  the 

CDMA  channel,  new  packet  transmissions  combine  with  packet  retransmissions  to  form  the 
composite  channel  traffic.  This  resulting  composite  channel  traffic  can  be  characterized  by 
its  intensity  and  steady-state  arrival  distribution.  The  composite  traffic  intensity  v  simply 
equals  the  sum  of  the  new  and  retransmitted  traffic  intensities  {u  =  t/o  +  t'r)*  The  composite 
arrival  distribution  /m(0  steady-state  probability  distribution  for  the  number  of 

attempted  transmissions  M  in  a  given  time  slot.  The  form  of  this  distribution  depends 
primarily  upon  the  relative  values  of  po  and  pr,  and  on  the  population  size  N.  Significant 
simplification  in  its  form  results  when  Po  =  Pr  =  P>  in  which  case  the  distribution  becomes 
binomial  with  parameters  p  and  TV.  Further  simplification  occurs  when  Po  —  Pr  0  nnd 
TV  — »  oo,  in  which  case  the  binomial  distribution  approaches  the  Poisson  distribution  with 
arrival  rate  equal  to  the  composite  traffic  intensity  v. 

A  Type  I  Hybrid  ARQ  protocol  is  used  for  error  control.  This  protocol  combines 
the  forward-error-detecting  (FED)  capability  of  plain  (conventional)  ARQ  with  the  use 
of  forward-error-correcting  (FEC)  codes  to  achieve  improved  network  performance.  An 
example  of  a  TVpe  I  Hybrid  ARQ  channel  (hereafter  referred  to  as  Hybrid  ARQ)  is  shown 
in  Fig.  2.  The  ARQ  portion  of  this  protocol  is  of  the  stop-and-wait  (SAW)  variety  where 
the  sender  waits  for  either  a  positive  (ACK)  or  a  negative  (NAK)  acknowledgement  of  the 
packet  just  sent  before  sending  another.  It  is  assumed  that  acknowledgements  are  made 
instantaneously  over  a  separate,  noiseless  return  channel.  With  this  assumption,  plain 
ARQ  is  equivalent  to  the  standard  slotted  narrow-band  ALOHA  channel  access  protocol.  If 
a  NAK  is  received,  then  the  sender  retransmits  the  current  packet  in  the  next  time  slot  with 
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probability  p,.  Other  ARQ  schemes  such  as  go-back-N  and  selective-repeat  can  be  used  to 
achieve  higher  throughput  efficiency  but  at  the  expense  of  increased  network  complexity. 
The  FEC  portion  of  the  |>rotocol  is  used  to  combat  the  effects  of  poor  channel  conditions 
and  tends  to  reduce  the  number  of  retransmissions. 

For  this  analysis,  we  assume  that  the  channel  conditions  (i.e.,  background  noise 
power  spectral  density  (PSD)  excluding  the  multiple-access  interference)  remain  constant. 
Otherwise,  the  Type  I  Hybrid  ARQ  protocol  becomes  less  efficient  as  the  channel  error 
rate  decreases  because  the  extra  error  correction  parity  bits  are  wasted.  For  time-varying 
channel  conditions,  it  is  better  to  use  the  Type  H  Hybrid  ARQ  protocol  which  uses  only 
FED  when  the  channel  conditions  are  good  (plain  ARQ)  and  uses  both  FED  and  FEC  when 
channel  conditions  are  poor  [5].  If  channel  conditions  degrade  beyond  the  capabilities  of 
Type  n  Hybrid  ARQ,  code  combining  techniques  (adaptive  Hybrid  ARQ)  can  be  employed 
to  successfully  recover  the  packet  [6,7]. 

Each  packet  that  is  transmitted  over  the  Hybrid  ARQ  channel  of  Fig.  2  is  comprised 
of  one  or  more  code  words.  Each  code  word  is  used  for  simultaneous  error  detection  and 
error  correction.  The  CDMA  channel  portion  of  Fig.  2  corresponds  to  the  coding  channel 
for  the  network.  For  hard  decision  decoding,  the  coding  channel  u  a  binary  symmetric 
channel  having  crossover  probability  Ps  and  with  the  following  cutoff  rate  and  capacity: 

Ro  =  1-logj  (n.y^4P5(l-P5)),  (2) 

C  =  l-H[Ps],  (3) 

where  H[P5]  is  the  binary  entropy  function  given  by 


H{Ps)  =  -Ps  logj  Ps-il-  Ps)  log2(l  -  Ps).  (4) 


Expressions  for  the  probability  of  symbol  error  Ps  can  be  obtained  by  manipulating  (2) 
and  (3).  For  the  cutoff  rate  case,  the  probability  of  symbol  error  is 


Ps  = 


(l->/r^) 

2 


(5) 


where 

a  =  2^“**  —  1, 


(6) 
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and  for  the  capacity  case 

P5  =  F-M1-C],  (7) 

where  is  the  inverse  operator  of  H{').  The  probability  of  symbol  error  Ps  depends 

upon  the  type  of  modulation  used  by  the  physical  links  in  our  system.  For  this  (malysis, 
we  consider  the  use  of  binary  modulation,  e.g.  BPSK  or  DPSK,  and  results  are  developed 
for  the  DPSK  case.  The  uncoded  bit  error  probability  in  the  presence  of  additive  white 
Gaussian  noise  (AWGN)  is  /\  =  f(Xt)  where  A(  is  the  bit  energy>to-total  noise  ratio.  For 
DPSK,  /(A«)  =  |exp{— At},  and  for  BPSK  /(Aj)  =  Q{y/2X^))  where 

When  coding  is  used,  the  above  code  symbol  error  probability  becomes 

Ps  =  /(rAt),  (9) 

where  r  is  the  code  rate. 

CDMA  is  a  spread  spectrum  multiplexing  technique  which  is  characterized  by  the 
use  of  high  rate  (many  chips  per  code  symbol)  pseudonoise  (PN)  spreading  sequences. 
Spreading  sequences  are  selected  from  code  families  having  good  (thumbtack)  autocorrela¬ 
tion  function  properties  for  accurate  acquisition,  low  cross-correlation  function  properties 
for  reduced  multiple-access  interference,  and  membership  to  a  large  family  size  to  support 
a  large  number  of  users.  It  is  assumed  that  the  period  of  the  spreading  sequence  is  much 
greater  than  the  code  symbol  duration.  It  is  also  assumed  that  the  spreading  sequences 
are  completely  random  so  that  the  resulting  multiple-access  interference  can  be  modeled 
as  a  Gaussian  random  variable.  In  practice,  spreading  sequences  are  pseudorandom,  but 
for  certain  code  families  (Gold  codes)  the  Gaussian  assumption  is  valid  to  within  a  scale 
factor  [8,9].  For  this  analysis,  a  receiver-based  spreading  sequence  distributicm  is  used.  To 
transmit  a  packet,  a  sender  modulates  each  packet  bit  with  the  spreading  sequence  of  the 
intended  receiver. 

The  following  simplifying  assumptions  are  made: 

1.  The  usual  bursty  user  assumption  (pr  >  Po)  i*  relaxed  to  allow  these  probabilities  to 
vary  freely  over  their  entire  range  of  values  (0  <  Pr,Po  <  !)• 
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2.  The  traffic  characteristics  are  homogeneous  throughout  the  network.  Therefore,  the 
probability  distribution  of  the  multiple-access  interference  is  the  same  at  each  receiver. 

3.  For  the  N  =  oo  cake,  the  composite  traffic  arrival  distribution  is  Poisson  with  an 
average  arrival  rate  of  v  packets  per  slot.  For  the  N  <  oo  case,  the  comp>osite  arrival 
distribution  is  derived  below  by  means  of  a  Markovian  model. 

4.  All  users  employ  the  same  error  detecting/correcting  code. 

5.  All  users  employ  a  common  carrier  frequency,  code  symbol  duration  Ts,  and  modu¬ 
lation  technique. 

6.  All  users  employ  the  same  processing  gain  q  =  Ts/Tc,  where  Tc  is  the  PN  chip 
duration. 

7.  Each  user  has  acquired  coarse  and  fine  synchronization. 

8.  The  effect  of  multiple-access  interference  consists  of  additional  broad-bauad  noise  as 
discussed  earlier.  Therefore,  the  equivalent  one-sided  noise  PSD  due  to  trans¬ 
mitter  at  the  receiver  matched  to  the  desired  signal  is  Nj  =  '^jPTc,  where  7y  is  the 
attenuation  of  the  signal  and  P  is  the  transmitted  power.  The  attenuation  can  be 
used  to  account  for  fading,  or  users  having  unequal  transmitted  power. 

9.  The  user’s  signals  are  mutually  noncoherent  so  that  their  PSD’s  add. 

10.  Because  the  channel  is  not  a  capture  channel,  all  packets  arrive  at -the  ireceiver  at  the 
same  time  and  at  equal  power  levels. 

Finally,  we  will  represent  the  total  additive  noise  at  the  front  end  of  the  receiver 

as  (9] 

Nt  =  No-^-fi~''EEsij,  (10) 

>=> 

where  No  is  the  background  noise  PSD,  Es  is  the  code  symbol  energy,  and  m  is  the  number 
of  simultaneous  users.  Since  Es  =  rEi,  the  bit  energy-to-total  noise  ratio  for  the  user 
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(11) 


j.  _  _ 7n^» _ 

‘  +  »?-*»•£?»  E7=»  7,- 

i^n 

_  _ 7n^o _ 

1  +  JZ%i  7y  ’ 

where  Ao  is  the  bit  energy>to-background  noise  ratio.  With  the  above  equal  power  assump¬ 
tion,  (11)  becomes 


1  +  ri-^rXo{m  -  1) ' 


III  Network  Analysis 


In  this  section,  the  steady-state  and  dynamic  performance  of  the  Type  I  Hybrid  ARQ  pro¬ 
tocol  in  slotted  direct-sequence  CDMA  systems  is  analyzed.  Expressions  for  the  asymptotic 
throughput  (hereafter  referred  to  as  the  throughput)  and  corresponding  average  packet  de¬ 
lay  are  developed  for  the  steady-state  performance.  To  evaluate  the  dynamic  performance, 
expressions  are  developed  for  the  input-output  packet  flow  rates  of  the  network. 

m-A  Steady-State  Throughput-Delay  Analysis 

The  throughput  T  at  the  output  of  the  Hybrid  ARQ  channel  is  defined  as  the  expected 
number  of  successful  packets  S  per  slot; 

T  =  E{S}.  (13) 

By  the  chain  rule  of  expectation,  (13)  can  be  expressed  as 


where 


T  =  E{E{S  I  M}}, 


E{S\M}  =  Yikps\M{k\l) 

*=i 


P5|M(fc  I  0  ==  Prob{5  =  k  successful  packets  |  Af  =  /  attempted  transmissions), 

=  (t) (0|i  -  Pc(OI'‘‘- 


(16) 


Here,  Pc(0  is  the  probability  of  a  correctly  received  packet  when  there  are  /  simultaneous 
packet  transmissions.  Using  (15)  and  (16),  the  throughput  expression  in  (14)  becomes 

T  =  Z  |E  *  (*)  ^c(01i  -  Pc(/)]‘-‘|  MO,  (17) 

where  /mO)  is  the  probability  distribution  function  for  the  number  of  attempted  transmis¬ 
sions  (composite  arrivals)  during  a  particular  slot  interval.  Since 


{i;  ‘  (1)  = iPc(‘), 


T  =  Z‘Pc(')M‘)-  (») 

1=1 

Among  other  factors,  throughput  depends  upon  the  channel’s  code  rate  r,  processing 
gain  t},  and  the  bit  energy-to-background  noise  ratio  Ao.  By  combining  equations  (9)  and 
(12),  the  following  expression  can  be  obtained; 


where  /“^(O  is  the  inverse  operator  of  /(•)  and  Ps  is  the  code  symbol  error  probability 
given  by  (5)  and  (7)  for  the  cutoff  rate  and  capacity  cases,  respectively.  As  in  [3],  assume 
that  the  code  used  has  the  property  that  the  packet  error  probability  Pe{1)  =  1  -  Pc(0  ^ 
zero  for  r  <  and  one  if  r  >  {H;o,C).  For  the  cutoff  rate  case,  this  actually  occurs 

when  a  very  long  constraint  length  convolutional  code  is  used  with  sequential  decoding. 
Similarly,  for  the  capacity  case,  this  occurs  as  the  coding  complexity  becomes  arbitrarUy 
large.  In  either  case,  it  is  assumed  that  the  coding  overhead  due  to  forward  error  detection 
(FED)  is  small  and  can  be  neglected.  For  the  limiting  case  as  r  -t  {iZo,C}  and  for  fixed 
values  of  N,  tj,  and  Ao,  (20)  gives  the  maximum  number  of  simultaneous  packet 

transmissions  as  m  =  [mj .  If  the  actual  number  of  packets  /  exceeds  m,  then  information 
is  being  transmitted  at  a  rate  above  the  cutoff  rate  or  capacity  of  the  channel.  As  a  result, 
all  packets  are  incorrectly  received  and  must  be  retransmitted  (  Pc(l)  =  0,/  >  th  ).  If 
the  number  of  packets  is  less  than  or  equal  to  m,  then  all  packets  are  received  successfully 
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(  Pc{l)  =  1,1  <  m  ).  The  throughput  expression  in  (19)  becomes 

inin(fH,N) 

T=  £  I/mO).  (21) 

isl 

In  order  to  compare  throughputs  of  different  systems  on  an  equal  basis,  we  need  to  account 
for  {/Zo,C}  and  q.  The  normalized  throughput  or  network  utilization  becomes 


T{pr,{Ro,C},t,) 


<=i 


(22) 


This  is  the  average  number  of  successful  packets  (per  slot)  per  unit  time  per  unit  bandwidth. 
Note  that  the  normalized  throughput  is  a  function  of  the  probability  of  retransmission,  the 
code  rate,  and  the  processing  gain,  and  is  always  less  than  one.  The  normalized  throughput 
is  also  a  function  of  the  number  of  users,  the  composite  traffic  intensity,  and  the  bit  energy- 
to-background  noise  ratio.  However,  this  dependency  is  not  shown  explicitly  because  these 
parameters  are  assumed  to  be  uncontrolled.  The  maximum  normalized  throughput  is 


max  max  T(pr,{fIp,C},f?), 

0<Pr<l  0<{Jia<C)<l 

i<n 


max  max 

0<pf<l  o«Jt«,o)<i 
»<•» 


{igo,C} 

V 


min(fft,Ar) 

E  '/«(')• 

J=1 


(23) 


For  the  cutoff  rate  case,  (23)  represents  the  practically  achievable  upper  limit  on  the  max* 
imum  normalized  throughput  Tr,.  For  the  capacity  case,  (23)  represents  its  theoretically 
achievable  (best)  upper  limit,  Tc-  The  optimal  probability  of  retransmission,  code  rate, 
and  processing  gain  are  denoted  by  p*,  {Rl,C*},  and  respectively. 

The  composite  packet  arrival  distribution  /m(0  i®  derived  by  using  a  Markov  model 
similar  to  that  developed  in  [10].  The  channel  w  viewed  as  a  discrete-time  system  where 
X(t)  represents  the  number  of  backlogged  users  n  during  a  particular  slot  interval.  The 
network  may  then  be  in  any  one  of  N  +  1  possible  states.  A  geometric  distribution  for  the 
probability  of  retransmission  results  in  the  retransmission  delay  also  being  geometricaUy 
distributed  with  an  average 

Dr^Ejpr,  (24) 
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where  E  is  the  average  number  of  packet  retransmissions  per  successfully  transmitted 
packet.  The  memory  less  property  of  the  geometric  distribution  allows  for  a  simple,  sin¬ 
gle  state  description.  Assuming  N  and  po  to  be  time-invariant,  X{t)  represents  a  Markov 
process  with  state  transition  matrix  P  =  [pi^-]  where  the  state  transition  probabilities  are 
defined  as 

Pij  =  Prob{X(t  +  1)  =  i  I  X{t)  =  i}.  (25) 

Calculation  of  the  above  state  transition  matrix  is  more  complicated  than  in  the  narrow- 
band  case  [10],  but  is  simpler  than  in  [1]  because  Pc{l)  is  binary  in  value  (Pc(0  =  ^ 

I  <  m  and  Pc{l)  =  0  fot  I  >  m  ).  From  this  simplification,  the  following  expressions  for  the 
one-step  state  transition  probabilities  can  be  derived  (App>endix  A); 

(for  j  <  i) 

0,  i  -  j  >  m 

^  _  l,p^)6(l  -  j,  l,p,),  i-j<m 

(for  j  =  i) 

Hk,  N  -  i,p,)b{0,i,pr),  i  <  m 

+  EU<h+i  Hk,i,Pr)K^,N  -  I,p„),  m<i<  N 
Z^^^lHk,N,Pr)  +  biO,N,Pr),  i=N 


Pij 


=  < 


(26) 


(for  J  >  i) 

0,  y  <  y  -  I  <  m 

6(y  -i,N  -  i,Po)[l  -  H*.*.Pr)]. 

[  b{j  -i,N  -  i,po),  m<j<NJ-i>m 

where  the  terms  in  the  summation  represent  the  binomial  distribution  function  such  that 


6(a,n,p)  =  (a)p®(l  -p)* 


(27) 
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When  m  =  1,  (26)  reduces  to  the  standard  narrow-band  slotted  ALOHA  case  as  in  [10]. 
For  the  infinite  population  model,  if  we  let  W  — »  oo  and  Po  — »  0  such  that  Npo  = 
then  (26)  becomes 


Pa 


(for  j  <  i) 

0,  i  -  j  >  m 

(forj  =  0 

+  EU*fi+i  -  i,Po),  m<i  <  N 
Z^=^+ih{KN,pr)  +  6(0.Ar.p,),  i  =  N 


(for  y  >  ») 


0, 


(28) 


r>-»  --r 


y  <  "*»  y  “  *  ^  »w 

m  <  y  <  iv,  y  —  I  <  fh 

(  '  ’  m  <  y  <  AT,  y  -  I  >  rh 

With  P,  we  can  solve  for  the  equilibrium  state  probability  distribution  ir  in  the  following 
equation: 

£  =  £P,  (29) 


where 


1=  [ir(0),jr(l),...,»(Ar)].  (30) 

Using  the  above  solution  for  £,  we  can  now  express  /Ar(0  terms  of  the  number  of  back- 
logged  users  as 


/a/(0  =  IT  Ml  I  n)x(n), 


(31) 


n=0 
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where 


imn(i,iv-n; 

fM{t  1  n)  =  ^  h{k,  N  -  n,po)b{l  -  k,  n,pr).  (32) 

fc=niax(i— n,0) 

Substituting  (32)  into  (23),  the  normalized  throughput  equation  becomes 

r  p  inin(ifc,N)  N  adn{l,S-n) 

=  E  'E>'W  E  l(t.Ar-«,P.)4((-*.n.P,). 

l<i|  t=i  •*=0  ifc=m*x(i-n,0) 

(33) 

For  the  special  case  when  N  =  oo,  po  =  pr  =  0  and  PoN  =  (33)  becomes 

MRo,cy  =  max  max  ~ — ■ — -  >  - - i  r34) 

o<p,<i  o<{«o.c>si  17  /! 

i<ij  »=i 

which  is  the  same  as  the  throughput  expression  derived  in  [3]  where  FEC  is  employed 
without  retransmissions. 

The  average  packet  delay  D  ot  a  network  is  often  modeled  as  the  sum  of  a  random 
delay  component  Dry  a  deterministic  component  D^y  and  a  one  slot  transmission  time: 

D  =  Dr  +  Di  +  l.  (35) 

The  deterministic  component  D^  can  be  used  to  model  any  fixed  delay  inherent  in  the 
network  such  as  the  round-trip  delay  time  required  by  each  acknowledgement  when  using 
the  stop-and-wait  ARQ  scheme.  For 'this  analysis,  it  is  assumed  that  the  deterministic 
component  is  zero.  As  such,  equation  (35)  becomes 

D  =  Dr  +  l.  (36) 

The  random  delay  component  D,  is  the  average  retransmission  delay  defined  earlier  and 
represents  the  average  rescheduling  delay  or  the  average  backlog  time  experienced  by  a 
user.  This  random  component  is  necessary  to  prevent  backlogged  users  from  contending 
with  other  users  blocked  during  the  same  time  slot.  FVom  Little’s  result  [ll].  Dr  can  be 
expressed  as 


Dr^ 
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(38) 


where  fi  is  the  average  number  of  backlogged  users  which  can  be  calculated  from 


ft  =  53  n»(n). 

n=0 

By  using  (36)  and  (37),  the  throughput  can  be  directly  related  to  the  average  packet  delay 
(35)  as 

T<«.o)  =  (39) 

This  alternate  expression  for  the  throughput  is  the  one  that  is  actually  used  in  the  maxi¬ 
mization  required  by  (33).  The  throughput  expression  in  (23)  provides  a  convenient  check 
during  the  analysis. 

Throughput  can  also  be  related  to  E,  the  average  number  of  packet  retransmissions 
per  successfully  transmitted  packet,  by  the  following  expression 

V 


E  = 


-1. 


(40) 


where  ^  average  number  of  times  a  packet  must  be  transmitted  until  it  is 

successfully  received.  This  expression  can  be  used  with  the  data  from  section  FV  (Fig  4)  to 
obtain  a  plot  of  E  against  Substituting  (40)  and  (37)  into  (24),  the  probability  of 

retransmission  can  be  expressed  as 


n 


(41) 


A  similar  expression  for  the  probability  of  original  packet  transmission  can  be  obtained  as 

Vo 


Po  = 


-  n)w{n) 
T 

N-li 


(42) 


These  expressions  for  p,  uxl  Po  used  in  computing  the  P  matrix  (26)  and  are  key  in 
performing  the  maximization  in  (23). 


m-B  Dynamic  Throughput>Delay  Analysis 

The  analysis  presented  so  far  is  not  sufficient  to  fully  characterize  the  performance  of  our 
CDMA  network.  Up  to  this  point,  we  have  assumed  that  the  network  is  in  equilibrium. 


IS 


That  is,  for  a  given  number  of  backlogged  users  n,  the  input  transmission  rate  Tin  i*  ^ual 
to  the  output  or  delivered  transmission  rate  Tout-  The  channel  input  transmission  rate  is 
defined  as  the  number  of  new  packet  transmissions  per  time  slot.  The  normalized  input 
transmission  rate  is 

Tin=^^^{N-n)po.  (43) 

n 

Note  that  if  we  take  the  expectation  over  n  in  (43),  then  T,>,  becomes  the  input  traffic 
intensity  Vo-  The  channel  output  transmission  rate  is  defined  as  the  probability  of  up  to  m 
successful  packet  transmissions  in  a  given  time  slot.  The  normalized  output  transmission 
rate  is  given  by 

Tout  =  -  nyPo)h{l,n,Pr)  -  (1  -  Po)^"'*(l  -  PrT-  (44) 

Note  that  if  we  take  the  expectation  over  n  in  (44),  then  the  output  transmission  rate 
becomes  the  normalized  throughput  in  (22).  When  N  and  po  are  stationary  with  time,  (43) 
and  (44)  are  referred  to  as  the  channel  load  line  and  equilibrium  contour,  respectively. 

A  channel  is  said  to  be  stable  when  its  load  line  intersects  (nontangentially)  the 
equilibrium  contour  in  exactly  one  place.  Points  where  the  load  line  intersects  the  equi> 
librium  contour  are  defined  as  channel  operating  points.  For  stable  channek  with  a  large 
population  size,  analysis  has  shown  that  networks  actually  operate  close  to  these  channel 
operating  points  [12].  In  Fig.  3,  (44)  is  plotted  against  three  dififerent  channel  load  lines 
(43)  to  show  the  effects  of  changing  network  operating  conditions  on  the  network’s  stability. 
Load  line  A  represents  a  stable  channel  having  an  operating  point  (7x,  ni)  with  relatively 
high  throughput  and  a  low  number  of  backlogged  iisers  (low  delay).  Load  line  C  also  repre¬ 
sents  a  stable  channel,  but  in  this  case  the  network  is  overloaded  wiih  a  channel  operating 
point  having  relatively  low  throughput  and  a  high  number  of  backlogged  users  (high  delay). 
Load  line  B  represents  an  unstable  channel  having  three  possible  channel  operating  points 
with  the  most  desirable  operating  point  being  at  (rs,ns).  Assuming  that  the  network  is 
currently  operating  in  the  vicinity  of  this  point,  it  will  not  remain  there  because  X{t)  is 
a  random  process.  That  is,  there  is  a  nonzero  probability  that  the  number  of  backbgged 
users  will  exceed  na.  When  this  occurs,  Tin  >  Tout  which  causes  the  operating  point  to  drift 
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(accelerate)  up  the  load  line  to  the  overloaded  channel  operating  point  (low  throughput  and 
high  delay).  While  there  is  also  a  nonzero  probability  that  the  channel  operating  point  will 
return  below  the  {T3,n2}  point,  simulations  show  that  the  channel  tends  to  remain  near 
the  overloaded  operating  point  instead  [12]. 

For  the  finite  population  case,  if  the  channel  is  stable,  then  the  network  operating 
point  can  reside  in  the  vicinity  of  the  channel  o|>erating  point  indefinitely.  For  an  un¬ 
stable  chwnel,  however,  network  operation  around  the  desired  channel  operating  point  is 
achievable  for  only  a  finite  amount  of  time  before  the  channel  becomes  overloaded.  Finite 
population  channels  can  be  stabilized  by  decreasing  pr  (at  the  expense  of  increasing  delay), 
by  decreasing  N,  or  by  employing  external  stabilizing  measures  which  include  various  re¬ 
transmission  control  strategies  [13].  For  the  infinite  population  case,  channels  are  inherently 
unstable  and  can  operate  near  their  desirable  operating  poitns  for  only  a  finite  amount  of 
time  [10,14].  This  operational  time  can  be  increased  by  thi-  methods  mentioned  above  to 
achieve  performance  close  to  the  theoretical  optimum  [13]. 

IV  Performance  Evaluation 

In  this  section,  the  dependency  of  the  network’s  throughput-delay  performance  on  N,  u,  A^, 
Pri  {fZa>  and  rt  is  examined  in  detail  for  both  the  steady-state  and  dynamic  cases.  The 
optimed  probabilities  of  retransmission,  processing  gains,  and  code  rates  are  identified.  The 
results  obtained  here  for  DS/DPSK  signaling  can  also  be  obtained  for  DS/BPSK  signaling 
within  the  same  framework. 

rV-A  steady- state  Throughput-Delay  Performance 

Figs.  4  and  5-6  are  plots  of  the  maximum  normalized  throughtputs,  Tr,  and  Tc,  against 
trafiic  intensity  v  and  against  average  delay  D,  respectively.  These  plots  illustrate  the 
combined  effects  of  maximizing  throughput  over  the  probability  of  retransmission  pr,  and 
the  code  rate  and  processing  gain  The  first  maximization  required  by  (23)  is 

performed  by  fixing  D  in  (39),  allowing  pr  to  vary  (0  <  pr  <  1)  in  (29),  and  solving  for  a 
(ft,  x(n))  solution  which  maximizes  (39).  This  procedure  is  repeated  for  all  possible  values 
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of  rh  (1  <  fh  <  A')  and  results  in  N  optimal  (Pr>7{jt^c})  pairs  for  each  fixed  value  of  D.  An 
alternate  procedure,  whereby  the  throughput  is  fixed  and  the  delay  is  minimized,  produces 
the  same  optimal  (pr *”,2))  combinations.  For  the  second  maximization  required  by 
(33),  D  is  again  fixed  for  a  given  level  of  background  noise.  The  code  rate  and  processing 
gain  {Ro,C-,tj}  are  allowed  to  vary  in  (20)  yielding  m  values  for  which  the  corresponding 
throughputs  are  known  from  step  one.  Each  value  of  throughput  is  then  normalized  and 
the  maximum  normalized  value  is  selected.  Note  that  (20)  requires  a  solution  for  Ps-  For 
the  cutoff  rate  case,  a  solution  for  Ps  is  found  in  a  strught  forward  manner  by  using  (5) 
and  (6).  For  the  capacity  case,  a  Ps  soluticm  is  found  by  solving  (3)  and  (4)  simultaneously. 

Results  in  Figs.  4  and  5-6  are  shown  for  population  sizes  N  =  10, 30,  and  oo;  and 
for  background  noise  levels  of  A<,.=  oo  and  Ao  =  10.34  dB.  The  background  noise  level  of 
10.34  dB  was  chosen  as  it  results  in  a  bit  error  probability  of  10~^  for  uncoded  DPSK  in 
the  absence  of  multiple-access  interference.  Although  the  infinite  population  case  may  be 
an  unrealistic  population  size,  it  does  provide  a  useful  bound  on  system  performance. 

IV'A.l  Optimization  of  Throughput-Delay  Performance  over  Probability  of 
Retransmiasion 

Optimal  values  for  probability  of  retransmission  exist  for  all  values  of  trafiic  intensity  as 
shown  in  Fig.  7.  For  the  low  range  of  trafiic  intensities  {v  <  1.0),  there  are  two  distinct  p* 
values  for  each  value  of  v.  Selecting  the  lower  p*  value  for  a  given  trafiic  intensity  results  in 
the  maximum  possible  value  for  {Tr.,Tc}  (Fig-  4)  but  with  nonminimal  delay  (Figs.  5-6). 
In  Figs.  5-6,  these  lower  p,  values  correspond  to  the  upper  branch  of  the  throughput-delay 
curve.  Here,  v  is  increasing  along  the  curve  from  A  to  B  and  back  to  C.  Points  on  the 
A-B  portion  of  the  curve  correspond  to  stable  operating  points,  whereas  those  on  the  B-C 
portion  of  the  curve  do  not.  A  maximal  {7'jt.,Tc}  value  is  achieved  under  the  following 
conditions:  p*  C  1.0,  p„  —*  1.0,  and  A  — »  A.  As  a  result,  the  throughput  is  closely 
approximated  by  the  equation  {Tr^,Tc}  =  v  for  low  i/  (i/  <  0.5)  in  Fig.  4,  and  by  the 
equation  =  N/D  for  the  upper  branch  (A-B)  in  Figs.  5-6.  These  two  throughput 

expressions  are  related  through  Little’s  result  (37)  by  applying  the  above  conditions  to  (41) 
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and  (39),  respectively. 

Selecting  the  higher  p*  value  for  a  given  traffic  intensity  {v  <  1.0)  in  Fig.  7  results  in 
a  nonmaximal  {Tjt„,Tc}  value  (not  shown  in  Pig.  4  )  but  at  minimal  delay  (Figs.  5-6).  The 
higher  values  correspond  to  the  lower  right-hand  branch  of  the  throughput-delay  curve 
in  Figs.  5-6  (increasing  v  from  E  to  F).  This  throughput-delay  curve  agrees  with  the  results 
from  earlier  throughput-delay  analyses  for  conventbnal  narrow-band  slotted  ALOHA  [12]. 
For  these  nonmaximal  {Tk^,Tc}  values,  p*  >  Po  and  6  is  relatively  low. 

Note  in  Fig.  7  that  for  the  low  range  of  traffic  intensities  (i/  <  1.0),  the  upper  and 
lower  curves  for  p,  are  separate  and  distinct.  For  a  given  traffic  intensity,  this  separation 
represents  the  change  in  p,  that  is  necessary  to  move  the  network  from  an  operating  point 
having  maximal  throughput  to  one  having  minimal  delay,  and  vice-versa.  In  Fig.  5,  for 
example,  the  throughput-delay  points  for  traffic  intensities  of  v  =  0.525  and  1.00  are  plotted 
for  comparison  purposes  for  the  two  values  of  p*.  Note  in  this  example,  that  for  i/  =  1.0 
it  seems  that  an  improvement  in  throughput  can  be  had  for  a  modest  increase  in  delay  by 
operating  on  the  upper  branch  (B-C).  However,  the  upper  branch  (B-C)  operating  point  is 
not  stable,  whereas,  the  lower  branch  (E-F)  operating  point  is.  Stability  issues  are  discussed 
further  in  section  IV-B.  Note  that  for  this  range  of  v,  p*  values  do  not  depend  on  A^. 

For  the  upper  range  of  traffic  intensities  {u  >  1.0)  in  Fig.  7,  there  are  also  two 
possible  optimal  values  for  the  probability  of  retransmission  p*  as  indicated  by  the  solid 
and  dashed  lines.  The  solid  line  corresponds  to  the  case  when  either  plain  or  Hybrid 
ARQ  is  the  optimal  protocol.  Here,  p*  can  be  approximated  by  p*  =  u/N.  Because  p*  is 
independent  of  A^,  the  same  value  of  p*  is  used  for  both  plain  and  Hybrid  ARQ  for  a  given 
traffic  intensity.  The  dashed  line  corresponds  to  the  case  when  Type  1  Hybrid  ARQ  with 
CDMA  (hereafter  referred  to  as  CDMA)  is  used.  Here,  p*  =  v/N  exactly  and  the  composite 
arrival  distribution  /a/(/)  is  binomial.  Note  in  Fig.  7  that  the  (dashed)  p*  =  u/N  lines  are 
shown  for  their  widest  possible  range  of  values.  Where  they  actually  begin  (lower  left  end) 
depends  on  the  cutoff  intensity  (discussed  below)  which,  in  turn,  depends  on  Ao.  Where  the 
dashed  lines  end  (upper  right  end)  depends  on  the  populatbn  sise  and  the  traffic  intensity. 
At  high  intensity,  high  processing  gain  causes  rh=  N.  As  a  result,  there  are  no  packets  to 
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retransmit  and  the  p*  =  v/N  lines  terminate. 

In  Fig.  4,  the  plain/Hybrid  ARQ  p*  values  correspond  to  the  plain/Hybrid  ARQ 
curves  and  the  CDMA  p*  values  correspond  to  the  upper  branches  which  break  away  from 
the  plain  ARQ  curves  at  various  triiffic  intensities.  In  Figs.  5-6,  the  B-C-A  branch  of 
the  throughput-delay  curve  (1/  increasing  from  B  to  C  then  toward  A)  corresponds  to  the 
plain/Hybrid  ARQ  p*  values  and  the  C-D  branch  corresponds  to  the  CDMA  p*  values. 

Note  in  Figs.  4  and  5-6  that  for  the  plain/Hybrid  ARQ  case,  throughput  becomes 
vanishingly  small  and  delay  becomes  unbounded  as  the  traffic  intensity  increases.  How¬ 
ever,  when  CDMA  is  used,  this  degradation  in  performance  does  not  occur.  Instead,  the 
maximum  normalized  throughput  (Fig.  4)  decreases  to  some  minimal  value,  at  a  particular 
traffic  intensity,  and  then  increases  again.  The  corresponding  delay  decreases  (Figs.  5-6). 
The  traffic  intensity  at  which  CDMA  achieves  this  improved  performance  over  plain/Hybrid 
ARQ  is  defined  as  the  cutoff  traffie  intenaity  and  is  denoted  as  {t/^,  1/^}.  In  Fig.  4,  the 
cutoff  intensities  occur  at  the  intersection  of  the  CDMA  curves  with  the  plain/Hybrid  ARQ 
curves.  Note  that  for  Ao  =  00  and  u  <  Tjt„  and  Tc  are  identically  equal.  For 

Ao  <  00  and  u  <  {*'**,1'®},  ^  Tc,  but  these  differences  are  usually  too  small  to  show 

graphically,  the  only  exception  being  for  Tr,  at  Ae  =  10.34  dB.  In  Figs.  5-6,  point  C  is  an 
example  of  a  throughput-delay  point  which  corresponds  to  a  cutoff  intensity.  Note  how  the 
delay  drops  precipitously  (C-D)  when  CDMA  is  used.  Recall  that  the  above  throughput- 
delay  results  must  be  adjusted  for  any  fixed  delay  inherent  in  the  particular  network  being 
considered.  Additional  fixed  delay  causes  the  curves  in  Figs.  5-6  to  be  shifted  upward  by  a 
corresponding  amount. 

IV-A.2  Optimization  of  Throughput-Delay  Performance  over  Code  Rate  and 
Processing  Gain 

An  important  parameter  arising  from  the  maximization  of  the  normalized  throughput  over 
the  processing  gain  and  code  rate  in  (33)  is  the  cutoff  traffic  intensity,  described 

as  follows.  1£  V  <  then  the  optimal  processing  gain  is  unity  for  all  bit 

energy-to-background  noise  ratios  (Ao  <  00).  If  i'  >  then  >  1  with 
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the  actual  value  depending  on  the  trafiBc  intensity,  population  size,  and  the  bit  energy-to- 
background  noise  ratio  as  described  belosr. 

Fig.  8  illustrates  the  effect  of  bit  eaergy>to-background  noise  ratio  A0  and  population 
size  N  on  the  cutoff  traffic  intensity  In  general,  {*/**>,  increases  to  its 

maximatn  value  (N)  as  Ao  decreases.  Note  that  for  each  there  is  a  corresp>onding 

Ao  which  is  defined  as  the  noise  limit  for  that  cutoff  intensity,  denoted  by  {Af”,Af}. 
Also,  for  each  Ao  there  is  a  corresponding  cutoff  traffic  intensity.  The  {6.82  dB,  4.81  dB} 
noise  liniits  are  defined  as  the  network  asymptotic  noise  limits  and  are  denoted  {A^,  Af }. 
The  network  cannot  operate  at  bit  energy-to-background  noise  ratios  below  the  asymptotic 
noise  limits,  because  they  represent  the  smallest  Ao  that  can  be  present  for  reliable  coded 
communications  in  the  absence  of  multiple-access  interference  [3]. 

Fig.  9  shows  how  the  optimal  processing  gain  depends  on  the  bit  energy- 

to-background  noise  ratio  Ao  for  different  traffic  intensities  u  and  population  sizes  N.  In 
general,  the  optimal  processing  gain  increases  with  decreasing  Ao.  For  the  jV  =  00  case, 
increases  toward  infinity  as  Ao  decreases  to  the  asymptotic  noise  limit  {A^,  Af }. 
For  the  iV  <  00  case,  decreasing  Ao  causes  a  similar  increase  in  but  toward  a 

finite  maximum  value  at  the  noise  limit  corresponding  to  a  given  traffic  intensity.  For  a 
particular  traffic  intensity,  if  {A*®,Af}  <  A*  <  {Af®,Af}  then  the  optimal  processing  gain 
is  unity.  Note  that  for  most  values  of  Ao,  depends  primarily  on  1/  and  not  so 

much  on  the  population  size. 

Figs.  10-11  show  how  the  optimal  processing  gain  depends  on  the  traffic 

intensity  for  the  different  population  sizes  N  and  bit  energy-to-background  noise  ratios  Ao. 
For  the  N  =  00  case,  {>7^,  *7^}  depends  on  v  in  a  nearly  linear  fashion.  The  effect  of  finite 
population  size  is  a  'staircase’  type  curve  which  tends  to  the  N  =  00  curve  as  N  — »  00. 

Fig.  12  summarizes  the  dependency  of  the  optimal  code  rate  {/2*,C*}  on  the  bit 
energy-to-background  noise  ratio  X„  and  population  size  N.  Note  that  the  selection  of 
the  optimal  code  rate  depends  primarily  upon  1/  and  Ao,  and  does  not  directly  depend  on 
N.  When  Ao  =  00  and  1/  <  the  optimal  code  rate  is  unity.  Hence, 

the  optimal  protocol  is  plain  ARQ.  In  fact,  plain  ARQ  is  used  only  when  Ao  =  00  and 
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V  <  When  A®  <  oo  and  v  <  {0.48,0.50}  <  (22*, C*}  <  1.0  In  this 

case,  Hybrid  ARQ  is  used.  For  Ao  <  oo  and  u  >  {i/*®,  v^},  the  optimal  code  rate  {22*, C*} 
is  {0.48,0.50},  respectively.  CDMA  is  the  optimal  protocol  in  this  case. 

It  should  be  noted  that  for  the  entire  range  of  traffic  intensities,  it  is  possible  to  oper¬ 
ate  CDMA  with  a  delay  D  =  1.0.  Unity  delay  is  achieved  by  selecting  the  appropriate  code 
rate  and  processing  gain  combination  so  that  m  =  iV.  In  most  cases,  however,  unit  delay 
is  achieved  at  the  expense  of  a  far  from  optimal  throughput.  In  general,  the  throughput 
at  unit  delay  TjOsi  is  very  small  for  0  <  »/  <  {i/®*,*/^}.  At  u  =  {*'*“>  J'^},  there  is  a  step 
increase  in  Td=x,  after  which  Td=\  — »  T^r„,c}  tav  N. 

Finally,  note  that  for  a  given  level  of  background  noise,  the  optimal  code  rate  can  be 
detemuned  as  the  traffic  intensity  varies  by  using  Figs.  8  and  12  .  For  a  given  Ag,  Fig.  8  gives 
the  corresponding  {i/®®,i/^}.  As  the  traffic  intensity  varies  above  and  below  this  particular 
{tA®<>,i/^},  the  optimal  code  rate  is  given  by  the  lower  branch  and  upper  branch  of  Fig.  12, 
respectively. 

IV-B  Dynamic  Throughput*Delay  Performance 

When  the  traffic  intensity  of  our  network  rises  to  a  level  at  or  above  the  cutoff  intensity, 
the  employment  of  CDMA  stabilizes  what  would  otherwise  be  an  unstable  network.  This 
effect  is  shown  in  Fig.  13  for  a  traffic  intensity  v  =  3.0,  a  population  size  of  N  =  10,  and  no 
background  noise  (A^  =  oo).  The  cutoff  intensity  in  thb  case  is  v®®  =  2.85.  Note  that  the 
plain  ARQ  curves  display  an  unstable  characteristic  with  an  overloaded  channel  operating 
point  K,  whereas  the  CDMA  curves  display  a  stable  characteristic  with  a  singular  channel 
operating  point  J.  Average  operating  points  generated  by  numerical  analysis  are  shown  for 
each  case.  Points  G  {Tr„  =  0.129,ft  =  7.417)  and  H  (Tij,  =  0.151, fi  =  2.70)  correspond 
to  the  plain  ARQ  and  CDMA  cases,  respectively.  For  the  CDMA  case,  note  that  the 
relatively  large  distance  between  the  average  operating  point  H  and  the  channel  operating 
point  J  can  be  attributed  to  the  small  population  size.  In  Fig.  5,  throughput-delay  points 
corresponding  to  i/  =  3.0  and  AT  =  10  for  both  the  plun  ARQ  (Tji^  =  0.129,  D  =  56.0)  and 
CDMA  {Tr^  =  0.151,  D  =  19.0)  cases  are  shown  for  reference.  In  Fig.  6,  a  throughput-delay 
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point  for  u  =  3.0  and  iV^  =  30  is  also  shown.  Observe  that  for  the  {Tjt„,Xo  =  oo)  curve, 
there  is  no  1/  =  3.0  throughput-delay  point  on  the  CDMA  curve  (C-D)  because  1/  =  3.0  is 
less  than  the  cutoff  intensity  for  N  =  30,  A<,  =  00  (v'*®  =  3.17).  Thus,  for  a  given  traffic 
intensity  and  Ao,  an  increase  in  N  causes  a  corresponding  increase  in  {1/^,1/^}  which,  in 
this  case,  destabilizes  the  network. 

Finally,  note  that  CDMA  achieves  a  similar  stabilizing  effect  for  Ao  <  00  and  v  > 
{t/^,1/^}.  For  this  case,  Hybrid  ARQ  operates  at  some  code  rate  less  than  imity  depending 
on  the  background  noise  level.  In  the  above  example  {N  =  10, 1/  =  3.0),  if  Ao  is  high  enough 
so  that  iZo  1.0  then  Ro  is  close  to  unity  and  is  unity  in  (43)  and  (44)  and  the  resulting 
Hybrid  ARQ  stability  curves  are  very  close  to  those  shown  for  the  plain  ARQ  case  in 
Fig.  13.  That  is  ,  the  Hybrid  ARQ  case  exhibits  the  same  unstable  characteristic  (ie., 
multiple  network  operating  points).  If  15.25dB  <  Ao  <  00  then  Ro  <  1.0  and  f)  =  1.0  and 
the  resulting  equilibrium  curves  generally  shift  to  the  left  but  maintain  the  same  unstable 
form.  For  the  CDMA  case,  15.25dB  <  Ao  <  00  causes  the  CDMA  stability  curves  to  shift 
in  a  similar  manner  while  maintaining  the  same  stable  form. 

V  Operational  Considerationg 

The  results  presented  above  suggest  that  the  network  will  operate  in  one  of  three  modes: 
plain  ARQ,  Hybrid  ARQ,  or  CDMA.  The  mode  of  operation  selected  depends  primarily 
upon  the  traffic  intensity  (/,  the  level  of  background  noise  Ao,  and,  to  a  lesser  degree,  the 
population  size  N.  All  of  these  quantities  must  be  monitored  by  each  user  for  proper 
network  operation.  For  the  case  when  the  traffic  intensity  is  low  enough  to  permit  stable 
network  operation  (  eg.,  v  <  0.5  ),  the  network  will  operate  in  the  plain  ARQ  mode  when 
Ao  =  00  and  in  Hybrid  ARQ  mode  when  Ao  <  00.  In  a  practical  system  with  finite  packet 
lengths,  the  plain  ARQ  mode  could  be  used  for  a  sufficiently  high  Ao.  In  this  case,  a  user  will 
have  to  decide  whether  it  wants  to  operate  with  optimal  delay  (upper  p*  values  in  Fig.  7) 
or  optimal  throughput  (lower  pi  values  in  Fig.  7)  by  choosing  the  appropriate  probability 
of  retransmission.  While  in  either  mode  of  operaticm,  the  probability  of  retransmission  is 
adjusted  according  to  variations  in  the  traffic  intensity  and  the  population  size  (Fig.  7)  and 
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the  code  rate  is  adjusted  according  to  variations  in  the  background  noise  level  (Fig.  12). 

When  the  traffic  intensity  reaches  a  certain  level  (0.5  <  u  <  1.0),  the  network 
becomes  unstable.  As  the  traffic  intensity  continues  to  increase,  a  cutoff  intensity  is  even¬ 
tually  reached  where  the  optimal  protocol  becomes  CDMA.  In  order  to  determine  when  to 
switch  to  CDMA,  an  upper  traffic  intensity  threshold  is  set  at  some  value  between  (/  ftt  0.5 
and  v  =  This  threshold  point  will  most  likely  be  set  closer  to  v  ts  0.5  th2m 

u  =  to  allow  the  system  time  to  change  over  to  the  CDMA  mode  of  operation. 

Note  that  a  lower  threshold  traffic  intensity  must  also  be  set  within  the  v  <  0.5  traffic 
intensity  range  to  trigger  the  return  to  the  plain/Hybrid  ARQ  mode.  During  the  CDMA 
mode  of  operation,  the  probability  of  retransmission  is  adjusted  according  to  variations  in 
the  traffic  intensity  and  population  size  (Fig.  7),  and  the  processing  gain  is  adjusted  ac¬ 
cording  to  changes  in  the  traffic  intensity,  the  background  noise  level,  and  the  population 
size  (Figs.  8-11).  The  code  rate  remains  constant  (Fig.  12). 

Finally,  it  should  be  noted  that  the  code  rate  and  processing  gain  will  be  chosen 
so  that  the  communication  link  will  be  operating  at  the  cutoff  rate  (capacity)  if  and  only 
if  the  actual  number  of  packets  in  a  slot,  I,  is  equal  to  m.  If  m  is  relatively  large,  then 
Prob(/  =  m)  may  be  small  so  that  the  system  operates  either  above  or  below  the  cutoff 
rate  (capacity)  most  of  the  time.  Hence,  the  asymptotic  performance  bounds,  especially 
those  in  terms  of  the  cutoff  rate,  may  be  ‘more  achievevable’  in  CDMA  systems  than  in 
narrow-band  systems  where  m  =  1. 

VI  Concluding  Remarks 

This  paper  has  presented  an  analysis  of  the  steady-state  and  dynamic  performance  of  Hy¬ 
brid  ARQ  protocols  in  slotted  direct-sequence  CDMA  channels.  Numerical  analyses  were 
presented  to  show  how  network  design  parameters  such  as  the  probability  of  retransmission, 
code  rate,  and  processing  gain  should  be  chosen  in  order  to  maximize  system  performance. 
It  was  also  shown  how  the  above  design  parameters  depend  on  the  uncontrolled  network 
parameters  such  as  population  size,  traffic  intensity,  and  bit  energy- to-background  noise  ra¬ 
tio.  Important  parameters  such  as  the  cutoff  intensities,  noise  limits  and  asymptotic  noise 
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limits  arising  from  this  analysis  helped  illustrate  when  it  is  best  to  employ  one  of  the  various 
error  control  techniques  discussed  in  this  paper  (ARQ,  Hybrid  ARQ  ,  or  CDMA).  Many 
important  properties  were  revealed  that  are,  quite  likely,  general  characteristics  of  CDMA 
networks.  These  include  the  general  nature  of  the  throughput  versus  traffic  intensity  and 
throughput  versus  delay  characteristics,  the  dependency  of  the  probability  of  retransmis¬ 
sion  on  the  traffic  intensity  and  the  critical  role  that  the  probability  of  retransmission  plays 
in  the  throughput-delay-stability  tradeoff,  the  linear  dependency  of  the  optimal  processing 
gain  on  the  traffic  intensity,  the  dependency  of  the  optimal  code  rate  on  the  bit  energy-to- 
background  noise  ratio,  the  relative  insensitivity  of  the  system’s  performance  to  changes  in 
the  number  of  transceivers,  and  the  stabilizing  effects  of  employing  CDMA. 
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A  Appendix 

In  this  appendix,  the  state  transition  probabilities  p,-y  are  derived  for  the  Markov  model 
discussed  in  section  III- A.  Recall  that  we  view  the  channel  as  a  discrete-time  system  where 
X{t)  represents  the  number  of  backlogged  users  n  during  a  particular  slot  interval.  For  a 
population  size  of  N  paired  transceivers,  the  network  may  be  in  any  one  of  N  +  1  possible 
states.  The  state  transition  probability  Pij  is  defined  as  the  probability  of  moving  from  a 
state  having  t  backlogged  users  in  a  given  time  slot  to  a  state  having  j  backlogged  users  in 
the  next  time  slot  or 

p.y  =  Prob{X(f  +  1)  =  i  I  X{t)  =  i}.  (45) 

In  one  step,  X{t)  may  decrease  by  an  amount  up  to  rh,  remain  the  saime,  or  increase  in  size 
up  to  N. 

Case  I  :  A  decrease  in  backlog. 

An  I  =  i  ~  j  decrease  in  occurs  when  /  retransmissions  are  successful.  In 
the  present  system,  this  can  occur  only  when  I  retransmissions  occur  together  with  up  to 
th  —  t  new  transmissions.  Any  combination  of  new  transmissions  greater  than  —  /)  plus 
/  retransmissions  would  be  blocked  by  the  channel.  Any  combinations  involving  fewer  or 
greater  them  /  retransmissions  would  represent  a  different  state  transition  probability.  Thus, 

Pi,i-i  =  Prob[up  to  (m  —  /)  NTX]  •  Prob(/  RTX] 

-  51  KKff  -i,Po)h{l,i,Pr),  (46) 

where  the  terms  in  the  summation  represent  the  binomial  distribution  function  such  that 

t(a,n,p)=(S)p«(l-p)"-,  (47) 

and  NTX  and  RTX  are  abbreviations  for  new  packet  transmisaions  and  packet  retrans¬ 
missions,  respectively.  Note  that  m  limits  the  size  in  which  decreases  in  the  number  of 
backlogged  users  can  be  made.  That  is 

Pi,i-l  =0,  for  /  >  m  (48) 
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Case  n  :  Backlog  remains  constant. 

Remaining  at  the  same  number  of  backlogged  users  can  occur  in  three  ways.  IVans- 
missions  may  involve  neither  NTXs  nor  RTXs,  only  NTX,  or  only  RTXs.  Note  that  trans¬ 
missions  involving  both  NTXs  and  RTXs  necessarily  increase  or  decrease  X(t). 

A.  A  transmission  involving  neither  NTXs  nor  RTXs  occurs  with  probability 

Pii  =  Prob[no  NTX]  •  Prob|no  RTX] 

=  6(0,N-i,p„)K0,t,Pr).  (49) 

B.  A  transmission  involving  only  NTXs  can  maintain  the  same  number  of  backlogged 
users  only  if  the  number  of  NTXs  are  less  than  or  equal  to  m.  This  occurs  with  probability 


Pii  =  Prob[no  RTX]  •  Probjup  to  m  NTX] 


inin(Ar— i.ifi) 

=  (1-Pr)*  E  Hk,N-i,p,).  (50) 

k=0 

C.  A  transmission  involving  only  RTXs  can  maintsdn  a  constant  backlog  size  if  more 
than  m  RTXs  are  attempted.  This  occurs  with  probability 


Pii  —  Prob[no  NTX]  •  Prob[more  than  m  RTX] 

=  b{0,N-i,po)  E  (51) 

Using  (49),  (50),  and  (51)  and  accounting  for  the  relative  sizes  of  the  current  number 
of  backlogged  users  t,  m,  and  iV,  p,-,-  can  be  expressed  as 


».»H)  ff  _  f,p^)J(o,i,p,.),  i<m 

+  Efc=«a+i5(*:,*,Pr)5(0,N-i,po),  m<i<N 


(52) 


Hk,  N,pr)  +  6(0,  N,pr),  i  =  N 
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Case  ni  :  An  increase  in  backlog. 

An  /  increase  in  the  number  of  backlogged  users  occurs  when  /  NTXs  occur  together 
with  more  than  (m  —  /)  RTXs.  This  occurs  with  probability 


Pii  =  Prob[/  NTX]  •  Prob[more  than  (if*  —  1)  RTX] 

=  (53) 

ksO 

Note  that  as  opposed  to  the  decrease  in  backlog  case,  a  one-step  increase  in  backlog  is  limited 
by  the  quantity  N  —  i.  Note  also  that  increases  in  backlog  are  impossible  for  combinations 
of  NTXs  and  RTXs  totaling  less  than  or  equal  to  th. 

Combining  (46),  (52),  and  (53),  the  state  transition  probabilities  pij  become: 

[  (for  3  <  0 


0,  i  —  3  >  m 

^  _  *,Po)5(i  -  i,t,Pr),  *  -  i  <  m 

(fori  =  *) 

^  _  *,p«)5(0,*,p,),  •  <  If* 

+  EUm+i  ^(*.*tPr)ft(0,  N  -  i,po),  m<i<  N 
bik,N,pr)  +  biO,N,pr),  i  =  N 


(54) 


(for  j  >  *) 


0,  J  <  »b,  j  -  I  <  m 

b{3  -i,N-  i,Po)(l  -  K*.*.Pr)),  m<3  <  N,3 -i  <m 

biJ-ifN -i,po),  m<3<N,3-i>m 


Note  that  for  in  s  1,  (54)  reduces  to  the  state  transition  probabilities  for  the  conventional 
narrow-band  slotted  ALOHA  [10],  and  that  (54)  is  a  special  case  of  the  state  transition 
probabilities  derived  in  [1]  for  random  access  CDMA  networks. 
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Figure  1:  CDMA  Network  Packet  Flow  Model 
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Figure  2:  Type  I  Hybrid  ARQ  Channri 
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Figure  3:  CDMA  Network  Stability  for  Different  Network  Operating  Conditions 
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Figure  4:  Maximum  Normalized  Throughput  versus  IVaffic  Intensity 
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Figure  5:  Average  Delay  versus  Maximum  Normalized  Throughput  for  N  =  10 
(N  =  30  results  are  overlayed) 
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OPTIMAL  PROBABILiry  OF  RETRANSMISSION,  p,* 


Figure  7:  Optimal  Probability  of  Retransmission  versus  IVaffic  Intensity 
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Figure  8:  Cutoff  lYaffic  Intensity  versus  Bit  Energy-to-Background  Noise  Ratio 
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Figure  9:  Optimal  Processing  Gain  versus  Bit  Energy-to-Background  Noise  Ratio 
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Figure  10:  Optimal  Processing  Gain  versus  l^affic  Intensity  for  the  Cutoff  Rate  Case 
(Capacity  Case  at  Ao  =  oo  is  shown  for  reference) 
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Figure  11:  Optimal  Processing  Gun  versus  IVafiic  Intensity  for  the  Capacity  Case 
(Cutoff  Rate  Case  at  Ao  ==  oo  is  shown  for  reference) 
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OPTIMAL  CODE  RATES 


Figure  12:  Optimal  Code  Rate  versus  Bit  Energy-to-Background  Noise  Ratio 
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Figure  13:  Number  of  Backlogged  Users  versus  Input/Output  Rate  Ekiuations  for  CDMA 
Network  Stability  Analysis 


